Role of electron correlations in transport through domain walls in magnetic nanowires 
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The transmission of correlated electrons through a domain wall in ferromagnetic quasi-one- 
dimensional systems is studied theoretically in the case when the domain wall width is comparable 
with the Fermi wavelength of the charge carriers. The wall gives rise to both potential and spin 
dependent scattering. Using a poor man's renormalization group approach, we obtain scaling equa- 
tions for the scattering amplitudes. For repulsive interactions, the wall is shown to reflect all incident 
electrons at the zero temperature fixed points. In one of the fixed points the wall additionally fiips 
the spin of all incident electrons, generating a finite spin current without associated charge current. 

PACS numbers; 73.63.Nm, 71.10.Pm, 75.70.Cn, 75.75.-|-a 



Introduction Electronic properties of magnetic wires 
with domain walls (DWs) attract much interest because 
of possible applications of the associated magnetoresis- 
tance effect So far, the puzzlingly huge magnetoresis- 
tance up to thousands of percents in Ni wires and struc- 
tured magnetic semiconductors is not well understood 
0, 0, 0] . Most of the existing theories of the DW resis- 
tance do not take into account the electron correlations 
1MB. 13. 13.1^ 1 except for the case of a wide DW (adiabatic 
regime) [lO|. However, in effectively one-dimensional 
(ID) structures one should take into account the ten- 
dency to the non- Fermi- liquid behavior possible en- 
hancement of the electron-impurity interaction |0| and 
of the localization effects jT^. Moreover, non-adiabatic 
DWs (smaller or comparable to the Fermi wavelength) 
can be achieved in ferromagnetic semiconductor wires 
with constrictions. 

In this Letter we study the effect of electron corre- 
lations on the resistance of a ferromagnetic wire with 
a non-adiabatic DW. Since the DW in a ID wire acts 
as a localized spin-dependent scattering center, a strong 
influence of electron correlations is expected [l^ . We 
apply the renormalization method [ij, [iSj for the loga- 
rithmically diverging terms in the perturbation theory of 
the electron interactions. We find the zero-temperature 
fixed points for repulsive interactions, in which the wall 
reflects all incident electrons. This might explain the gi- 
ant magnetoresistences observed in magnetic nanowires. 
Moreover, we show that for one fixed point the electron 
spin is reversed at the reflection. This leads to a nonzero 
spin current with no charge current, which is of high in- 



*On leave from Departamento de Fi'sica, Universidade de Evora, 
P-7000-671, Evora, Portugal 



terest for applications in spintronics devices. 

Model and method We consider a magnetic ID system 
with a local exchange coupling between conduction elec- 
trons and a spatially varying magnetization M(r). The 
wire itself defines the easy axis (2- axis), and a DW cen- 
tered at z = separates two regions with opposite mag- 
netizations: Mz{z ±00) = ±A/o. Assuming M(r) to 
lie in the xz plane and the DW width to be smaller than 
the Fermi wavelength, one can write the single-particle 
Hamiltonian as (J > 0): 

Ho = - — -r^ + hV5{z) + JAUz)a, + h\6{z)&^, (I) 
2m az^ 

where the term fi\5{z)(Tx describes spin dependent 
scattering due to the Ma:{z) component 0, ?iA = 
J Mx{z)dz, and V is a potential (spin independent) 
scattering term. We assume that spin-t electrons are spin 
majority ones for z < and spin-minority for z > 0. 

An incident electron from the left with momentum k 
and spin ] (or [) can be transmitted to the z > Q region 
preserving its spin, but changing momentum from k to 

k- (or fc+), given by k^ = ^k'^±AJMoni/n^. If the 
transmission occurs with spin reversal, the momentum k 
remains unchanged. The reflection amplitudes for spin-cr 
electrons with or without spin reversal are denoted by 
and r„, respectively. The same convention applies to 
the transmission amplitudes t'^ and t„. The scattering 
amplitudes are given by: 

^ta)(^) - (. + .T + 2^y)2 + 4A2 =-ffl)^^) + ^'^^) 

AiXv 

4(i)W = = ^Ta)('=)' (3) 

with V — Tik/m, = hk^/m, where the upper (lower) 
sign refers to t (!)• We shall henceforth denote by 



2 



e(±p, a) the energy of a scattering state with momentum 
+p (or —p) and spin a, incident from the left (or right). 
The scattering amphtudes satisfy general relations that 
can be found from a generalization of the Wronskian the- 
orem 0,^3 to spinor wave- functions: the Wronskian of 
any two states having the same energy is a constant, 

W{^Pi,^2) = ^l(^)^-^ M^) = ^onst, (4) 

where '0* denotes the transpose of the spinor ^. 

In order to deal with the electron interactions, it is con- 
venient to rewrite the scattering states in second quan- 
tization form, making use of right (dgo-) and left (bqa-) 
moving plane wave states. We introduce operators Ck,a 
and dk.a for the eigenstates corresponding to electrons 
incident from the left and right, respectively. The plane 
wave operators are linear combinations of the operators 
of scattering states. Electron interactions are then intro- 
duced through the Hamiltonian 

dkidq , 
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where the Greek letters denote spin indices, and the sum- 
mation convention over repeated indices is used. The 
coupling constants gi and (72 describe back and forward 
scattering processes between electrons moving in oppo- 
site directions, respectively. Because the Fermi momen- 
tum is spin dependent, we distinguish between gi(2)], 
which describe interaction between spin-majority parti- 
cles (that is spin-f on the left and spin-| on the right 
of the barrier) and 5i(2)Xj which describe interaction be- 
tween spin-minority particles. We use gn2)± to denote 
interaction between particles with opposite spin. The 
forward scattering process between particles which move 
in the same direction will not affect the transmission am- 
plitudes, although it will renormalize the Fermi velocity 
|llL H^ . This effect is equivalent to an effective mass 
renormalization and the electrons with different spin ori- 
entations may turn out to have different effective masses, 
in which case our calculations remain valid |l7l |. 

Following Ref . 0| , the corrections to the transmission 
amplitudes are calculated to first order in the perturba- 
tion Hint- Such corrections diverge logarithmically near 
the Fermi level and will be dealt with in a poor man's 
renormalization approach. The perturbative correction 
to ta{p') can be obtained from the perturbation expan- 
sion for the Matsubara propagator 

GAr) - ^{T^e~I"-^^^'^'^'d,Ar)cl,Jo, (6) 

where (...)o denotes the average in the Fermi sea of scat- 
tering states. The zero-order propagator for a =t is: 



g\ '(iLu) = 



1 



e(p',T) 



ztt(p') 



p — p' + iO P ~ p' 



where denotes a positive infinitesimal. The poles in the 
denominators identify the semi-axis on which the electron 
behaves as a right moving plane wave. The transmis- 
sion amplitude appears associated with the denominator 
p — p'~ — iQ which, for the variable p, gives a pole in 
the upper half plane. The meaning of this pole is that 
the transmitted particle is right-moving in the z > 
half-axis. Our strategy is to calculate the first order cor- 
rection term to in which a pole in p — p'^ — «0 will 
appear with the residue —ibt^(j)'\ which is the transmis- 
sion amplitude correction. 

The diagrammatic representation of q'^^^ is shown in 
Fig. The horizontal lines represent the electron scat- 
tered by the Hartree-Fock potential of the Fermi sea. 
Only processes where the electron is back-scattered by 
the Fermi sea are logarithmically large 0|. Consider, for 
instance, the upper left diagram - an electron, initially in 
state Cp',1 close to the Fermi level passes through the DW 
as a right-moving (a) particle. Then, it is reflected (from 
a to 6 particle) while exchanging momentum q with the 
Fermi sea on the z > semi-axis. Finally, it is reflected 
by the DW again, becoming a spin-up right moving par- 
ticle of momentum p. A logarithmic divergence occurs if 
the polarized Fermi sea can provide exactly the momen- 
tum that is required to keep the electron always near the 
Fermi level in the intermediate virtual steps. According 
to the physical interpretation of the diagrams, we always 
know on which side of the DW the interaction with the 
Fermi sea (closed loop in the diagram) is taking place. 
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FIG. 1: Feynman diagrams for the first order contribution 
C/'^' to the propagator @. The scattering state is represented 
by a double line, the a (b) particle is represented by a con- 
tinuous (dashed) line. The loop represents the Hartree-Fock 
potential of the Fermi sea. The scattered electron exchanges 
momentum q with the Fermi sea. 

We use Fermi level velocities v± and wavevectors kF± 
for majority or minority spin particles, henceforth. It can 
be seen that gi± terms are proportional to log|fci?_|_ — 
kp^\, so they do not diverge. Logarithmic divergence 
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would be restored in a spin degenerate system {kp+ = 
kp-). This can be understood from the diagrams in 
Fig. n as follows: the electron with spin a is reflected 
by a polarized Fermi sea with spin —a. The momentum 
provided by the Fermi sea is 2kF-on while the momen- 
tum required by the electron is 2kFa- The 172^ terms 
produce logarithmic divergences that would not exist in 
the absence of spin- flip scattering {t' — r' = 0). For the 
calculation of St'^ [p') the propagators we need to consider 
are -{Trap-a{T)c^^,^^). 

The expression for 5t^{p') is directly proportional to 
log(|e'|/D), where e' denotes the energy of the scat- 
tered electron and D is an energy scale near the Fermi 
level where the electron dispersion can be linearized 
[5t'\{p')/\og{\e'\/D) is given by the right-hand side of 
equation JHl below]. The logarithmically divergent per- 
turbation can be dealt with using a renormalization pro- 
cedure 101; reducing step by step the bandwidth D and 
removing states near the band edge is compensated by 
renormalization of t| . Applying this procedure and not- 
ing that <| + (5t| remains invariant as D is reduced, one 
finds the following differential equation: 

Now, we introduce the variable ^ = log(D/r'o), which is 
integrated from to logQe']/ Dq), corresponding to the 
fact that the bandwidth is progressively reduced from 
Do to |e'| (which will be taken as temperature: |e'| = T). 
It is convenient to rewrite the interaction parameters as 
gni) = (52T(i) -5iT(i))/4/iw+(-), g± = 92±/2h{v+ + V-). 
The scaling differential equations for the transmission 
amplitudes are 

^ = 5i [ rlrfy + rlr\t[ ] + [ r*r^t^ + r*r\t\ ] 



= '^girlr'fy + 25Tr|rT4 + 2g^ [r[*r^t^ + r\*r[t\ ] .(9) 

Equations for the reflection amplitudes r„{p') 
and r'^{p') can be obtained from the propagators 
-(T'rfop,±o-(T)Cp,^). The equation for r|(p') is 

+g^[r\*r[r^ +r'^r\r^ +r'{t^t\ ^r'^m\]-g^r^ (10) 
and the equation for r'| [p') is 

+9A. K*nn + + A*t\t\ + r'*t^t^] - g^r\ . (11) 



The scaling equations for spin-| amplitudes follow from 
the above by simply inverting the spin and velocity in- 
dices. Equations l|S|l- lfTT)) are the central result of this 
paper. Theorem Q gives V-/v+ = = ryr'y A 

standard second-order renormalization group treatment 
shows that in a ID magnetized system the coupling con- 
stants in equation ((S)) are not renormalized because the 
logarithmically divergent contributions cancel each other. 

Fixed points We have performed a numerical analysis 
of the scaling equations using the DW model O, with 
1^ = 0, for the initial parameters. We now analyze the 
nature of the fixed points predicted by the scaling equa- 
tions. The parameters of the model which enter the scal- 
ing equations are g^, g^, g±, and the ratio V-/v+. For 
repulsive interactions {g-\,gi,g± > 0) the system flows to 
insulator fixed points. For A/w+ larger than about 0.1, 
all transmission amplitudes vanish faster than any reflec- 
tion amplitude as T — > 0. In this limit we may rewrite 
the scaling equations for r-f, and the theorem ne- 
glecting the small transmission amplitudes. Theorem Q 
for W{^l^,ipk~,i) tells us that r|r^ + r^r[* = 0. The 
charge conservation condition is satisfied solely by the 
reflections, 
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from which we conclude that \r^\ = |?'x I the fixed point. 
Equations l|10|) and (|ll|l now take the form 

dr^' 
~d^ 
dr' 



(25^^gT-.9i)(l-kTl') (13) 



-g,-2g^)[l-—\r\\^)r',. (14) 



In the derivation of p3|l and 1)14(1 only the smallness of 
the transmissions amplitudes was assumed. The phases 
of the complex numbers r-f , are unchanged during scal- 
ing. The two fixed points we may consider correspond to 
r| approaching 0, or \r-^\ approaching 1. 

The situation \r^ \ — > requires 2g± — g] — gi > and, 
by charge conservation we have — > a/ v+ jv- . Upon 
integrating ((Tljl with ^ going from to log(r/I?o), the 
amplitude will vary from r'^ q to (T) . Introducing 
the reflection coefficient 'R!^^ = (w_/ti+)|r'| we obtain 



2(ST+9i-2s_L) 



(15) 



If 25r± - 5T - 3i > tl^cn VJ^ (T) ^ 1 as T ^ 0. The 
DW reflects all incident electrons while additionally re- 
versing their spin - it is a 100% "spin- flip reflector" at 
zero temperature, generating a finite net spin current but 
no charge current. In order to find the low T behavior of 
transmissions we put = 0, \r'^ \ ~ ^Jv^JvZ in the scal- 
ing equations for the transmission amplitudes and obtain 
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In the regime where gi+gi — 2g± > we have TZ'^{T) 
0, 7e|(r) 1. The DW reflects then ah incident elec- 
trons while preserving their spin. The scaling equa- 
tions for the transmission amplitudes yield \ ^ T^t+S-i ^ 

If 5T + 5i - 2gj. = then both r'^{T) and r„{T) tend 
to finite values. The scaling equations for t-^, t'^, with 
constant reflection amplitudes, become a linear algebraic 
3 by 3 system. The eigenvalues of the matrix give the 
temperature exponents and each transmission amplitude 
will be a linear combination of the three powers of T. For 
decreasing temperature there may be crossovers from one 
exponent to the other and the lowest one dominates as 
T ^ 0. 

For smaller values of A/u-|- {i.e., smaller than about 
0.1) in the Hamiltonian the system flows to a fixed 
point where r'^ vanishes about as fast as the transmissions 
and \ra \ 1. The transmission amplitudes still scale to 
zero as ~ TSt+Si, \t'^\ ^ T^f". The scaling equation 



for r' can be linearized in 



Neglecting the second 



order terms in t', and considering that |rcr| — > 1, we 
obtain |7^'^| - T'^iai+ai), in the scaling equation for 
r| we cannot neglect the terms containing transmission 
amplitudes on the right hand side. The behavior of Te- 
as T —> can be found from the charge conservation 

condition: 1 - \r^\^ ~ ymin{2(gT+Si):4g„}^ 

We may estimate the A parameter in Eq.(Q by as- 
suming that M(z) — Mocos9{z)z + Mosm9{z)x with 
cos9{z) = tanh(z/L), where L is the length of the DW 
[lo|. We then have A = irMoL/h, implying that 



2^ 



■Km _ n „ — (-l^kF^) 
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rfc 
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The condition for the DW to be smaller than the Fermi 
wavelength is Lkp+ < 2tt. The smaller Fermi wavelength 
is that of spin-majority electrons, 2t: /kp^. The ratio 
V-/v^ depends on the degree of polarization of the elec- 
tron system. In a ID nonmagnetic system there is a single 
Fermi momentum, hp, for up and down electrons and a 
Fermi energy Ep — h^kp/{2m). Once the system be- 
comes magnetized, there is a Zeeman energy shift of the 
bands, AE/2 = JMq, and the two new Fermi momenta, 
kp±, satisfy 



kp± AE 

= 1 ± 



kp AEp 
Inserting this result in Ea. H16() above, we obtain 



A _ (AE/AEp) 
v+ ^""[l + iAE/AEp)]' 



{Lk+) 



(17) 



(18) 



In the full polarization limit kp^ — 0, fcF+ ~ 2kp, 
and Eq. (fTH)l gives X/v+ w 0.79 Lkp+. Typical values 



for a non- fully polarized system are Ep = 90 meV and 
AE — 30 meV JJ. In this case we have v-/v+ = 0.84 and 
Eq. JIHl gives X/v+ K, 0.22Lkp+. Therefore, if Lkp+ is 
smaller than about 27r, the system can fiow to any of the 
fixed points described above. 

Lateral quantization may produce several channels. 
The possibility of inter-channel scattering then arises due 
to two causes: (i) electron interactions (which would 
require a modification of our theory to allow for inter- 
channel scattering); (ii) impurity scattering. For the lat- 
ter to be negligible the electron mean free path must be 
larger than the size of the constriction pinning the DW. 

The above spin-flip reflector DW was not found in Ref 
lOl . This is because the adiabatic DW considered in Ref 
inl | is a poor spin-flip reflector at the noninteracting level 
already - as in the regime of small \/vj^ above. 
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